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Abstract. Consider a frictionless surface S in a gravitational field that need 
not be uniform. Given two points A and B on S, what curve is traced out 
by a particle that starts at A and reaches B in the shortest time? This paper 
considers this problem on simple surfaces such as surfaces of revolution and 
solves the problem two ways: First, we use conservation of mechanical energy 
and the Euler-Lagange equation; second, we use geometrical optics and the 
eikonal equation. We conclude with a discussion of the relativistic effects at 
relativistic velocities. 



1. Introduction 

In 1696, Johann Bernouli posed the following "Brachistochrone problem" Find 
the shape of the curve down which a bead sliding from rest and accelerated by gravity 
will fall from one point to another in the least time. This problem assumes that 
the particle is falling on a vertical plane in a uniform gravitational field. Sir Isaac 
Newton, Gottfried Leibniz, Guillaume De L'Hospital, Jakob Bernouli, and Johann 
Bernouli showed that the solution is a cycloid, the curve traced out by a point on 
the rim of a rolling circle, see Dunham [lj. Solutions of the classical Brachistochrone 
problem typically use techniques of calculus of variations, see Gelfand and Fomin 
[2] , or geometrical optics, see Erlichson [3J . 

The problem of finding Brachistochrone curves with coulomb friction lying on a 
vertical plane in a uniform gravitational field has been discussed by Ashby et al. 
[4], Hayen [5], Heijden and Diepstraten [6], and further generalized to a curve with 
friction lying on a cylinder by Covic and Veskovic [7J. Vratanar and Saje [S] discuss 
the related problem of finding the brachistochrone curve in a non-conservative resis- 
tance field. The problem of finding Brachistochrone curves on cylinders in uniform 
fields has been solved by Yamani and Mulhem [9] and on cylinders and spheres and 
in the unpublished work of Palmieri 10 . The generalization to non-uniform fields 
has been discussed in the works of Aravind [11], Denman [12], and Venezian [15] . 
who find Brachistochrone solutions in linear radial fields, and the works of Denman 
[12] . Parnovsky [14], and Tee [15], who find solutions in inverse square radial fields. 
Further generalizations of this problem to include special relativistic effects have 
been studied by Farina [16] . Goldstein and Bender [17], and Scarpello and Ritelli 

This paper considers the following generalized problem: Let S be a smooth fric- 
tionless surface in a (not necessarily uniform) gravitational field. Given two points 
A and B on S, find the curve from A to B along which a particle released from A 
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reaches B in minimal time. In section 2 of this paper we use the Euler-Lagrange 
equation to present a theorem that solves the problem for a large class of surfaces in 
various gravitational fields. Our solution is general enough so that Yamani and Mul- 
herti' [9] solution on a cylinder is a direct application of our theorem as is Parnovsky 
[14j and Tee's |15j solution to the inverse square problem. Furthermore, generaliz- 
ing the result of Parnovsky [14] and Tee [15] , we prove that for a central force field 
with V = — r~ n no solution curve enters the sector — 2ir/ (n + 2) < 9 < 2ir/ (n + 2). 
In section 3 we generalize the optico-mechanical analogy used by Aravind [TTJ , Fa- 
rina [TB], and Parnovsky [T5] to surfaces. In particular, we use the eikonal equation 
[2T] to provide another proof of the theorem proved in section 2. Using the eikonal 
equation we find a solution for a charged particle moving at relativistic velocities 
lying in a uniform electric field which coincides with Goldstein and Bender's [17] 
and Farina's [16j result. Finally, we apply the light ray curvature equation |21j 
to prove that torsion always vanishes along Brachistochrone solutions in a central 
force fields. This result proves the implicit or explicit assumption taken by Aravind 
[TT] , Denman [12] , Venezian [13] , Parnovsky [14] , and Tee [15] that Brachistochrone 
solutions in radial fields are planar curves. 

The numerical computations and graphics in this paper were carried out and 
created in Mathematica [25], 

2. Solutions using Calculus of Variations 

Let S be a smooth frictionless surface in some gravitational field and let V (P) be 
the gravitational potential at point P on S. Choose distinct points A and B with 
V(A) > V(B) and assume that the velocity of a particle falling along a smooth curve 
from A to B is much less than the speed of light. Then by Newtonian conservation 
of mechanical energy we have 

1 fds\ 2 

o hr +v = v(A), 



2\dt / 

where ds is the arclength on S. Solving for (4f)^ gives 



2(V(A)-V)=(% 



2 

^ _> 

dt 



and separating variables yields 



dt = ± ^2(V(l)-V) dS - 



Thus the total time T is given by 



(2- 1 ) T = ± \i^i > ■ 




a \V(A)-V 

Now, let U be an open subset of M 3 and let x : U — > S be an orthogonal 
coordinate patch. Then x(u,v) = (xi(u,v),X2(u,v),X3(u,v)) is a subset of S in 
which ds 2 = Edu 2 + Gdv 2 , where E and G are the metric coefficents E = I s • 1^ 

1 au ou 

and G = • f^, and we may rewrite equation (|2.1[) in the form 
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and either as 

IT r B I e (—) 2 + c IT r B 

(2.3) ^ = ±V 2 A V -Yf^V- ^ = ^2 L ^ ^ ^ 

or 



IT f B IE + G ( —) 2 IT f B 

(2.4) T = ±] J- j A ^__i*L du = ± ^ ^ ^(u, 

where ""' denotes differentiation with respect to the other variable. 
To minimize T in (|2.3p . we must solve the Euler-Lagrange equation 

(2 5) ^ ^ ^ 

Since ^ = §g = = implies = (see 0), equation $23]) simplifies to 

— =C, 
<9it' 

where C is a constant that depends upon _B. By differentiating J- with respect to 
u' we obtain the separable differential equation 

= C. 

^{V{A)-V)(E{u') 2 + G) 
Squaring both sides and simplifying gives 



du _ / C 2 (V(A) - V)G 



dv I/ E(E-C 2 {V{A)-V))) 



whose solution expresses u as a function of v (see item 1 in Theorem [T] below) . 
Similarly, since ^ = when = 7^7 = ^7 = 0, we may use equation (|2 .4(1 and 
a modified form of (12.511 to obtain v as a function of u. 



Theorem 1. Let x : U — ► 5 &e an orthogonal coordinate patch on a smooth fric- 
tionless surface S. 

(1) J/ ^ = ^ = ^ = 7 then the solution to the Brachistochrone problem 
on S is given by the curve yi(u(v),v), where 



f v C 2 G(w)(V(A)-V{w)) 

(2) If ^ = ^ = ^ = 7 i/ien the solution to the Brachistochrone problem 
on S is given by the curve x(u,v(u)), where 



f U C 2 E(w)(V(A)-V(w)) 
2.1. Applications of Theorem [H 
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2.1.1. The Classical Brachistochrone Problem. Let us apply Theorem [T] to solve 
the classical Brachistochrone Problem. In this case, the particle is falling on the 
vertical plane given by x(tt, v) — (u, 0, v). Since E = G = 1 and F = 0, we obtain 




where V(v) = v and a is the particle's initial u-coordinate. Figure [1] illustrates 
several solution curves for a particle starting at the origin, each of which is uniquely 
determined by C and a choice of sign. Physically, the value of C is determined by 
the destination point B and the sign indicates motion to the left or right. Note 
that each of these curves represents a family of solution curves since the terminal 
point B can be positioned anywhere along the curve. 




Figure 1. Cycloid solutions in a uniform field 



2.1.2. Surfaces of Revolution. Let us consider the Brachistochrone problem on 
smooth surfaces of revolution S parametrized by the coordinate patches 

x n (tt, v) = (h(u) cosv, h{u) sin v, g(u)), (u, v) G Ix ((2n — 1)tt, (2n + 1)tt) 

y n (u,v) = {h(u) cosv, h(u) sinw, g(u)), (u, v) G Ix (2nir, (2n + 2)ir) , 

where g, h : R — > M are differentiable functions with h > 0, n G Z, and 

x„ (u, v) , n£ (2n7r, (2n + 1)tt) 
x n+1 (u,v), v G ((2n + l)7r, (2n + 2)7r) . 



Yn (U, V) 



Corollary 1. If S is in a uniform gravitational field parallel to the z-axis, the 
solution to the Brachistochrone Problem on S for a particle with initial position 
A = (h(uo), 0, g(uo)) is the curve 



7(V) 

where 



x„ (u, v (u)) , v (u) G [2nn, (2n + l)7r) 

y„ (u, v (u)) , v (it) G \{2n + 1)tt, (2n + 2)tt) , 




The right circular cone is obtained as a surface of revolution S by setting h (u) = 
g(u) = u > 0. For a particle with initial position A = (1, 0, 1), the solution curves 
7 (u) on S 1 given by Corollary [T] are defined in terms of 



v(u) = ± I 



I 2C 2 (uq-w) 
w 2 (w 2 — C 2 (uq — w)) ' 



where the sign ± positions the curve in one of the two half-spaces determined by 
the yz-plane. Figure [2] illustrates several solution curves. 
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Another interesting surface of revolution is the hyperboloid of one sheet obtained 
by setting h(u) — coshu and g (u) = sinhw. For a particle with initial position 
A = (1, 0, 0), the solution curves 7 (it) given by Corollary Q] are defined in terms of 

C 2 (sinhuo — sinhuj)((coshu>) 2 + (sinhw) 2 )) 
(cosh u>) 2 ( (cosh w) 2 — C 2 (sinhuo — sinhu;)) 

But unlike solution curves on the right circular cone, some solution curves on the 
hyperboloid do not attain a minimum z-coordinate but continue downward, spiral- 
ing around the hyperboloid and intersecting other solution curves as they do. When 
two spiraling solution curves intersect, one must evaluate the integral in (|2.2p to 
determine which curve minimizes time. 



v (it) = ± 





Figure 3. Solution curves on the hyperboloid of one- sheet 



2.1.3. Central Force Fields. Consider a particle confined to a plane and falling in a 
central force field with V = — \? for some n > 0. Using polar coordinates, we can 
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parametrize R 2 by x(u,v) — (ucos?;,Msinu), (u,v) £ (0,+oo) x (— 7r,7r); then the 
metric coefficients are simply E — 1 ,F = and G — u 2 . Now apply Theorem 1 
and set D = 1/C 2 to obtain the solution curves 



where A = x(mo, 0) is the initial position of A. Figure [4] illustrates several solution 
curves for an inverse square field with n = 1. In this case, the ± sign positions A 
in the first or fourth quadrant. 



Interestingly, when — ir < v < 0, numerical plots of solution curves asymptoti- 
cally approach the line 6 = ~2ir/3 as C — > 0. Since solution curves in the region 
< v < 7r are mirror images of those in the region — ir < v < 0, these calculations 
suggest that no solution curve enters the sector — 2it/3 < 9 < 2n/3. Note that 
solution curves seem to approach one of two limit curves that follow the rc-axis to 
the origin then continue along the rays 6 — —2tt/3 or 8 = 2ir/3. This is indeed the 
case and was observed and proved by Parnovsky [14] and Tee [15]. The following 
is an extension of their work for all n > 0. 

Without loss of generality assume that A = x(l,0) and the sign in (|2.8j) is 
positive so that the particle's trajectory lies in the fourth quadrant. Note that 
when 



(2.8) 






FIGURE 4. Solution curves for the inverse square field (n = 1). 
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the particle begins to move away from the origin and the sign in (|2.8| is reversed. 
We will show that equation (|2.9|) holds for exactly one positive real value by inves- 
tigating the zeros of the polynomial 

(2.10) r(w) = w n+2 D + w n - 1. 

Since r(0) = —1 and r(l) = D, the Intermediate Value Theorem tells us that 
r(co) = for some cq G (—1,1). Furthermore, 

r'(iu) = {n + 2)w n+1 D + nw™- 1 = w^ 1 ({n + 2)w 2 D + n) . 

But r'(w) > for all < w < D so r is increasing and cq is unique. Let 



(2.11) 



(1 - w n ) 



•j 2 (w n+2 D - (1 - w n )) 



dw. 



Since solution curves are symmetric in the line through x(co , 9) and the origin, we 
can determine the least upper bound of the polar angle 9 by evaluating its limit as 
cq — » 0. Using equation (|2. 10p we can express D in terms of Co by 



n — r ~ n ~ 2 - r~ 2 
u — c c • 



Then, equation (|2.11[) can be rewritten as 
(2.12) 9 = 



(1 - w n ) 



■ dw 



Co J 



1 — U) r 



w n + 2 (l ~ cl) - c" +2 (l - w n ) 
Let x = (f^)^; then (f2~T2")l reduces to 
2 M 



Cn 

n + 2 



1 



L + 2 

T - X i 



1 - CqX™+2 



\ c% +2 x-i(l - c») - c >+ 2 (l - c^) 



(fx 



n + 2 



n + 2 



1 - c x™+ 2 



o V \ (1 - c l ) - x 2 ( 



1 - cjjx ™+ 2 



(fa. 



l-c" 



■ dx 



\ ™+ 2 



Now, for c " < a; < 1, we have that 



l — cV:x n +' 2 



> 1. Therefore, 



n + 2 



1 -x 2 



■ (fa. 



c o \ mpr - * 

\ l-cffx ™ + 2 

Finally, it follows by an application of the Dominated Convergence Theorem that 



lim 9 = , 

c„^o n+2 J Vl - a; 2 



2 7T 

<fa = —7 (arcsin(l) — arcsin(O)) = 



n + 2 



(n + 2)' 
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Now, by symmetry the maximum polar angle of the particle is less than twice this 
limiting value. Therefore the limiting polar angle of a particle falling on a solution 
curve is -^ro and the central angle of the "forbidden sector" is 2it =rW = ^rk. 

n+2 ° n+2 n+2 

In summary, we have proved: 

Theorem 2. Let A{r,<p) E M 2 , r > 0. If V = — ^, n > 0, every Brachistochrone 
curve initiating at A lies in the sector <fr — < 9 < <j) + • 



3. Solutions using Geometrical Optics 

In this section we approach the Brachistochrone problem from the perspective of 
geometrical optics by considering the path of a light ray in a medium of nonuniform 
index of refraction. Let S : x (it, v) be a smooth surface in such a medium with 
index of refraction n (it, v). According to Fermat's Principle, light propagated from 
source A minimizes the total time T of travel to destination B, i.e., 

1 f B 

(3.1) T = — / n(u,v)ds, 

C . A 



where c is the speed of light in a vacuum. Note that when n(u, v) — y v(A)-v(u v) 
and c = lm/sec, equation (|3.1|l is equivalent to equation (|2.1|) . Therefore, com- 
puting the path of a light ray in this medium simultaneously produces a Brachis- 
tochrone solution. 

3.1. The Eikonal Equation. One can describe light rays in terms of their wave- 
fronts, which are level surfaces of some differentiable function L : M 3 — > M. More 
precisely, given a light ray a passing through a medium with index of refraction n, 
let T be the unit tangent. Then L is determined by the eikonal equation 

(3.2) V£ = nT. 



We can use (|3.2[) to solve the Brachistochrone Problem without appealing to the 
Euler-Lagrange equation. Assume that a is constrained to some smooth surface S 
with metric ds 2 = Edu 2 + Gdv 2 in some medium with index of refraction n(u,v). 
Suppose = f£ = 9n = and that 

(3.3) L = Cv + f(u). 

Then equation (|3 . 2|) yields the differential equation 

VL,VL = G- 1 C 2 +E- 1 [ =n 2 . 



du 



Solving this separable equation for / gives 




Using equations (|3.2p and (|3.3p . and letting e u = j^j and e^ = tt^-it we have 
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Expressing T in differential form gives 



idu ^idv \ G-G 2 /n 2 ^ , 
E2—e u + Gi— e v = J —4. — e u + CG ?e v 



ds ds J y G/r 

Therefore 

dvGi _ I C 2 /n 2 
d^E^ ~ Y G-C 2 /n 2 ' 
By solving this equation we obtain: 

Theorem 3. Let x : U — > S be an orthogonal coordinate patch on a smooth fric- 
tionless surface S. 

(1) If ^ = ^ = ^ = 0, i/ie light rays on S are given by x(u(v), v), where 




GC 2 In 2 

" (r) = ± I \l E(E-C 2 /n 2 ) dw - 
(2) If = ^ = ^ = 0, £/ie rays on 5 are given 6t/ x(u,v(u)), where 




< u ) =± l \l G{ G C -C 2 ln 2 ) ^ 

Corollary 2. Let x : U — > S 6e an orthogonal coordinate patch on a smooth 
frictionless surface S in a medium with index of refraction n 2 — V J_ V ■ 

(1) If Or- = 4^ = §^ = 0, we recover Theroem 1.1. 
v ' J au au au 7 

(2) // = = $f = 0, we recover Theorem 1.2. 

K / J av av av ' 

3.2. Special Relativistic Solutions. The strength of Theorem [3] lies in its abil- 
ity to give Brachistochrone solutions in situations where classical conservation of 
mechanical energy does not hold. For example, if the particle's velocity is near 
the speed of light, to minimize the time of travel in the lab frame (see [19]) we 
must replace the Newtonian mechanical energy equation by its special relativistic 
counterpart 

(3.5) 1C 2 ~c 2 + V = Vo, 

where c is the speed of light in a vacuum, V is the gravitational potential, and 
7 = (1 - (^f ) 2 c - 2 )-| (see [19]). Solving (33]) for 7 and then for (|f ) 2 gives 

'ds\ 2 c 2 (V -V)(V -V + 2c 2 ) 



dtj (Vo-V + c 2 ) 2 

From equation (13. ip it follows that 

fo c\ fl\ 2 _ (V -V)(V -V + 2c 2 ) 

[6 - b) [n] ~ (Vo-V + c 2 ) 2 ■ 

So, in sufficiently nice geometrical settings with reasonable potentials, we can apply 
Theorem [3] to find relativistic Brachistochrone solutions. 

For example, let us again consider a particle falling in a uniform gravitational 
field but confined to the vertical plane x(it, v) — (u,0,v). If the particle is falling 
at relativistic speeds we have 

1\ 2 = (vo-v)(vo-v + 2c 2 ) 
n J (vq — v + c 2 ) 2 
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where vq is the initial height of the particle. Note, in actuality the gravitational po- 
tential energy depends upon the relativistic mass and not the rest mass. Therefore, 
the potential we have chosen corresponds to a charged particle under the influ- 
ence of a uniform electric field (see Goldstein [19] and Goldstein and Bender [17]). 
Setting vq — and applying Theorem [3] gives the solution curves x(u(u), v), where 



u(v) 




k 2 /r 



1 - k 2 /r 



■ dw 



k 2 (2c 2 — w)w 



■ dw 



( A + 2c 2 (k 2 - l)w-{k 2 - l)w- 

and to avoid confusion we have indicated the constant of integration by k. Figure[5] 
illustrates several relativistic solution curves and classical Brachistochrone solutions 
with c = lOm/sec. The relativistic solutions are plotted in solid blue; the classical 
solutions are plotted with dashed green lines and closely approximate relativistic 
solutions in this region. However, Figure [6] illustrates situations in which such 
approximations are poor. 




Figure 5. Relativistic and classical Brachistochrone solutions. 




Figure 6. Relativistic and classical Brachistochrone solutions. 



3.3. Curvature of Light Rays. Finally, we consider the curvature of light rays in 
a medium with index of refraction n. Recall the following definitions and theorem 
from the classical theory of curves (see [20]h 

Definition 1. Let a(s) be a light ray with arc length parameter s and let T (s) be 

the unit tangent vector field along a. The curvature ofa(s) is defined by n — 

If n ^ 0, the unit normal and binormal vector fields along a are given by N = -^T' 

and B = T x N, respectively . 

Theorem 4. (Frenet-Serret) . If a(s) is a curve with arc length parameter s 
such that k(s) ^ for all s, then 

T' = kN, N' = -kT(s) + tB and B' = -rN. 

Definition 2. The function r(s) is called the torsion of a(s). 
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Now, light rays are distinguished curves with the following property: 

(3.7) Vn = nnN+ (T» Vn)T. 

In fact, when falling in a central force field, the torsion vanishes. To see this, 
differentiate (|3.7[) and get 

N' = (— \ nnN+ — [(Vn)' - k(T« Vn)N- ( T« Vn)'T] 
V kh J ku 



1 Y 1 
— ) Kn (Tm Vn) 

Kn I n 



N - — (T • Vn)'T + — (Vn)' . 
Kn Kn 



By Theorem |H the torsion of a light ray is determined by the component of (Vn)' 
in the direction of B. But for a particle falling in a central force field, n(r) = 
\Jr^r n I (t-q — r n ), where r is the distance from the origin and ro is the initial 
distance from the origin. Let e r denote the unit radial vector. Then, Vn = ^e r 

and (Vn)' = (4^) e r . Thus, (Vn)' and Vn are parallel and lie in the osculating 
plane. Consequently r = and Brachistochrone solutions in central force fields 
are planar. Futhermore, since Vn is parallel to the line from the particle to the 
origin, the osculating plane passes through the origin and contains the initial and 
terminal point of the curve. But, this is exactly the problem we investigated when 
we confined the particle to a plane. Therefore, a solution curve with critical radius 
Co and initial point A generates a surface of revolution whose meridians are solution 
curves. These remarks are summarized in our concluding theorem: 

Theorem 5. If a(s) is a solution curve lying in a central force field with critical 
radius Co, then r is identically and a(s) is a meridian of the surface of revolution 
generated by the curve W.< 
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